Thermodynamic and magnetic properties of the layered triangular magnet NaNiC>2 
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We report muon-spin rotation, heat capacity, magnetization, and ac magnetic susceptibility mea- 
surements of the magnetic properties of the layered spin-1/2 antiferromagnet NaNiC>2. These show 
the onset of long-range magnetic order below Tn = 19.5 K. Rapid muon depolarization, persisting 
from Tn to about 5K above Tn, is consistent with the presence of short-range magnetic order. The 
temperature and frequency dependence of the ac susceptibility suggests that magnetic clusters per- 
sist above 25 K and that their volume fraction decreases with increasing temperature. A frequency 
dependent peak in the ac magnetic susceptibility at T B f = 3 K is observed, consistent with a slowing 
of spin fluctuations at this temperature. A partial magnetic phase diagram is deduced. 

PACS numbers: 76.75.-K, 75.50.Ee, 75.30.Gw 
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Geometrically frustrated transition metal oxides ex- 
hibit a rich variety of magnetic behavior due to compet- 
ing interactions. A significant degeneracy in the ground 
state leads to such phenomena as spin glass 1 , spin liquid 2 , 
and spin ice^ phases. Triangular lattice antiferromagnets 
exhibit a variety of these phenomenal. When the trian- 
gles forming the lattice are distorted from equilateral to 
isosceles there is a partial release of the geometrical frus- 
tration, which can lead to more unusual forms of mag- 
netic order 5,6,7 . Among triangular lattice antiferromag- 
nets, LiNiO^, AgNiO^, and NaNi02 have somewhat 
enigmatic magnetic behavior. The difficulty of produc- 
ing stoichiometric LiNiC>2 has led to a variety of sample- 
dependent results^. AgNiC>2 can be produced in stoi- 
chiometric form but no magnetic Bragg peaks have so 
far been reported^. Recently neutron powder diffraction 
studies have determined the low temperature magnetic 
structure of NaNiO^ii*^- Darie et al— i find the ordering 
of the magnetic moments at 4 K to be a slight modifica- 
tion of the A-type antiferromagnetic ordering previously 
proposed**. The magnetic moments were found to be 
aligned at an angle of 100(2)° to the a-axis in the ac 
plane with no moment along the fe-axis. A peak in the 
magnetic susceptibility interpreted as the Neel temper- 
ature, Tn, has been observed around 20 KiiiiiS,. The 
Curie- Weiss constant, 6>cw = +36 K 15 , shows the pres- 
ence of ferromagnetic interactions above Tn- 

The intralayer and interlayer exchange constants of 
NaNi02, J\\ = — 13. 3 K and J± = 1.3 K, were determined 
from a model assuming an A-type antiferromagnetic or- 
dering with an anisotropy fieldi 6 -; the layers are suffi- 
ciently strongly coupled to permit long range magnetic 
order below Tn. The Ni-O-Ni bond angles are ~ 95° 
at room temperatureit An undistorted 90° geometry 
favours weak ferromagnetic superexchange, while a large 
deviation from a 90° bond angle can reverse the sign of 
this exchange coupling 1 —. In NaNiC>2 it appears that de- 
spite the distortion, in-plane ferromagnetic coupling pre- 
vails, though the precise nature of the spin and orbital 
ordering remains under discussioni2i22iSi. 

Above 480 K the space group of NaNi02 is rhombohe- 



dral (R3m) and there is a cooperative Jahn- Teller transi- 
tion to a low temperature monoclinic (C2/m) phase be- 
low this temperature 17 . The low-temperature structure 
can be considered to be layers of NiOe octahedra in the 
ab plane, with a trigonal distortion lengthening the Ni-0 
bonds in the ac plane along an axis at 41° to the c-axis 
of the crystal. The Ni 3+ (3d 7 ) ion is in the low spin state 
(t^gGg, S = 1/2), so the ground state is a singly-occupied 
|3z 2 — r 2 ) orbital with z along the axis of the Jahn- Teller 
induced trigonal distortion^. 

In this paper we present the results of zero-field muon 
spin rotation (//SR) experiments 22 on polycrystalline 
NaNiC>2 together with heat capacity and magnetic sus- 
ceptibility measurements. These provide evidence of long 
range magnetic order below Tn, a slowing of spin fluc- 
tuations at low temperature, and short range magnetic 
order within a small temperature region above Tn, with 
magnetic correlations persisting to higher temperature. 

NaNiC>2 was prepared from Na2C>2 and NiO powders 
heated at 700°C for 100 hours under pure oxygen flow, 
with intermediate grinding. X-ray powder diffraction 
showed that the impurity concentration was below the 
2 % resolution limit of the apparatus. Heat capacity data 
measured in magnetic fields between and 14 T, taken 
with a Quantum Design PPMS, are shown in Fig.^a). 
In zero field the transition at Tn = 19.5 K is seen as a 
rather broad peak, and there is no evidence for other 
phase transitions below 30 K. With increasing magnetic 
field the temperature of this peak decreases to about 14 K 
(see Fig. da)). 

Magnetic susceptibility data, taken using a Quan- 
tum Design MPMS SQUID magnetometer are shown in 
Fig. njb), (d) and (e). The high temperature dc sus- 
ceptibility data, shown in Fig. ^e), are consistent with 
$cw = +36 K 1 -. The real and imaginary parts of the 
ac susceptibility, %' and x", are presented in Figs, ^b) 
and (d) (driving field 3.5 Oe). They show a frequency de- 
pendent peak in %' located slightly above a spin-freezing 
temperature T s f , which is determined below. Figure [Jd) 
shows the temperature variation of this peak in \' with 
frequency. It was found that this could be fitted to the 
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FIG. 1: (Color online) The panels correspond to: (a) Heat 
capacity divided by temperature in fields between and 14 T. 
(b) Real part, x', of the ac magnetic susceptibility, (c) Tem- 
perature dependence of the peak in x' associated with Tf with 
a fit to the Ogielski relation (Eq. Q). (d) Imaginary part, 
x" , of the ac magnetic susceptibility, (e) Inverse of magnetic 
susceptibility data against temperature with a linear fit show- 
ing the high-temperature Curie- Weiss behavior. The vertical 
dashed lines indicate temperatures referred to in the text. 



Ogielski scaling relation^: 

T sf (l + (/ro) 1 /^), 



T : 



(1) 



where T is the temperature of the peak in x', / is the 
measurement frequency, To is the relaxation time of the 
system, z is a dynamic exponent, and v c is a critical ex- 
ponent. Fitting the data to Eq. Q gives T si = 3±0.2K, 
t = 5.4(2) x 10~ 3 s and zv c = 8.1 ± 0.4. This is typi- 
cal of glassy behavior or the slowing of spin fluctuations. 
Given that muon precession was observed down to 1.6 K 
(see below) , the formation of a true spin glass around T s f 
is excluded. One possible interpretation of the feature 
around T s f is that it could arise from a small concen- 
tration of defects comprised of pairs of Ni 2+ impurity 
spins which are associated with oxygen vacancies. Alter- 
natively, the fluctuations in the spins could be slowing 
down around T s f. At Tn, x' nas a small maximum and 
rises to a larger peak near 25 K, with the frequency de- 
pendence increasing from Tn to the peak, and decreasing 
above it. The temperature of the peak in x' &t 25 K 
decreases slowly with increasing frequency, which may 
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FIG. 2: (Color online) (a) Heat capacity divided by tempera- 
ture vs field at four temperatures, (b) Partial magnetic phase 
diagram deduced from heat capacity (•) and magnetization 
(o) data. AF: A-type antiferromagnetic phase. PM: Param- 
agnetic phase. FM: Ferromagnetic phase. 



be related to the presence of two sets of relaxation times 
varying differently with temperature (see below), x" a l so 
rises sharply near 25 K but to a plateau continuing up to 
~ 35 K. Together these suggest that short-range order 
persists up to ~ 25 K, that slowly fluctuating clusters of 
spins are present within a fast fluctuating paramagnetic 
bulk above this temperature, and that the volume frac- 
tion of clusters decreases with increasing temperature. 

Constant temperature heat capacity data with vary- 
ing magnetic field are presented in Fig.[3{a). Except for 
the data taken at 22 K, a peak is observed which cor- 
responds to the field labelled Hqx in the magnetization 
data reported in Ref. This suggests that this marks 
the upper field boundary of A-type antiferromagnetic or- 
der. At 22 K the heat capacity decreases with increasing 
field consistent with short-range order. The partial mag- 
netic phase diagram deduced from our heat capacity and 
magnetization data is shown in Fig.|2fb). 

Our zero-field //SR experiments were carried out us- 
ing the DOLLY instrument at the Paul Scherrer Insti- 
tute (PSI), Villigen, Switzerland. In our /iSR experi- 
ments, spin polarized positive muons (/i + , mean lifetime 
2.2 (txs, momentum 28 MeV/c) were implanted into poly- 
crystalline NaNi02 . The decay positron asymmetry func- 
tion, A(t) 22 , is proportional to the average spin polariza- 
tion of the muons stopped within the sample. The muon 
spin precesses around an internal magnetic field, P M , at 
a frequency — (7 p /27r)|P M |, where 7^/2^ = 135.5 
MHz T _1 . 

The asymmetry data were fitted to Eq. J2J 24 below Tn, 
and to Eq. © above Tn: 

A{t) = A(0) (T ie - Alt + P 2 e- A2i cos (27ri/ M t + fo)) (2) 
A{t) = A(0) (P { e~ Xit + P s e- X ^) , (3) 

where A(0) is the initial asymmetry. Pi and P2 are re- 
spectively the longitudinal and transverse components of 
the muon polarization, and P\ + P2 = 1. The exponen- 
tial relaxation associated with Pi reflects the dynamical 
fluctuations of the fields being probed. The P2 term de- 
scribes muon precession with a distribution of local fields 
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FIG. 3: Muon decay asymmetry in NaNi02 plotted at dif- 
ferent temperatures. The solid lines are fits of the data to 
Eqs. @ & © with the parameters shown in Fig. ^] 



ation is modelled with the two exponential components 
of Eq. |J3J , with the amplitude of the faster relaxing com- 
ponent decreasing with increasing temperature. Above 
24 K the relaxation is well described by a single expo- 
nential, Pfexp(— Aft), consistent with fast fluctuations of 
paramagnetic moments characterized by a single correla- 
tion time in the muon time window. 
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dephasing the muon spins. In a fully magnetically or- 
dered polycrystalline sample we expect P2/P1 = 2. Co- 
herent muon precession will be observed if long range 
order is present within the sample. P s and Pf describe 
slow and fast dynamic fluctuations respectively. A small 
initial phase offset, <fio, was observed below 2j\j, larger 
than could be attributed to errors in determining the 
time that the muons enter the sample. This could be 
produced by a small magnetic inequivalency in the po- 
sition of muons stopped within the sample, consistent 
with an asymmetric peak seen in Fourier transforms of 
the data. In the fitting procedure, data were fitted in 
the time range < t < 8 /is, where the effect of back- 
ground counts could be reliably subtracted. Rapid dy- 
namic fluctuations lead to Ai cx 7 M 2 (AP) 2 /^, where AB 
is the amplitude of the fluctuating local field and v is the 
fluctuation rate^i. 

Spectra measured at four temperatures are shown in 
Fig. [3J There are three distinct temperature regions ap- 
parent from the muon asymmetry spectra. At low tem- 
peratures (T < 19.5 K) there are clear oscillations in the 
asymmetry showing that long range magnetic order ex- 
ists and the observed ratio of Pi : P\ « 2 (see Fig. 0Jb)) 
indicates that the sample is magnetic over its entire vol- 
ume. The value of A2, s is much larger than Ai f (see 
Fig. ^c) & (d)) so at short times only the effect of A2, s 
is seen in Fig. [3J An intermediate temperature range 
(19.5 < T < 24 K) gives no oscillations, and the relax- 



FIG. 4: (Color online) Temperature dependence of the pa- 
rameters determined from fitting data to Eqs. (|5J & ©: (a) 
the oscillation frequency, u^, and the internal magnetic field, 
B^, with a fit to Eq. 0. (b) Amplitudes of the relaxation 
components Pi and P2, and Pf and P s . (c) Relaxation rates A2 
and A s . (d) Relaxation rates Ai and Af. The vertical dashed 
lines indicate temperatures referred to in the text. 

The temperature dependence of the parameters de- 
rived from fitting Eqs. J3J) & @ to muon asymmetry 
spectra are presented in Fig. The muon precession 
frequency, v tll in the ordered phase is shown in Fig.^fa). 
This is proportional to the sublattice magnetization at 
the muon site, and was fitted to a function 25 : 



u M {T) = ^(0)(1 - (T/T N )Y 



(4) 



The fit gives f M (0) = 64.2(2) MHz corresponding to a 
field at the muon site of ~ 0.5 T. Dipole field calculations 
show that this field will be experienced by muons near 
any of the oxygen atoms in the octahedron surrounding 
a nickel atom, in regions of high electron density 2 ^, and 
show that our results are consistent with the magnetic 
structure determined by Darie et al^i. Our calculations 
also suggest that the muon precession frequency is insen- 
sitive to small deviations from this magnetic structure. 
Fitting Eq. Q to the muon precession frequencies gave 
T N = 19.51(1) K and m = 0.24(1). This value of m sug- 
gests that the system is behaving as a 2D XY magnefe2I. 
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In relation to the peak just above T s f in the mag- 
netic susceptibility, we note that extrapolating the Ogiel- 
ski scaling relation to the muon time window suggests a 
maximum in the dynamic relaxation rate Ai should be 
observed around 7K, and a broad maximum of low am- 
plitude is just detectable at this temperature (Fig.QJd)). 
The presence of two exponential relaxation components 
above Tn (see Eq. suggests that short-range mag- 
netic order persists over a small temperature range of 
~ 5 K above Tn . The fluctuations in the magnetic field 
producing the faster relaxing component, A s , are two or- 
ders of magnitude slower than those producing the slowly 
relaxing component, Af. P s decreases with increasing 
temperature up to 24 K, showing that the ratio of slow to 
fast dynamic relaxation is decreasing. Above this tem- 
perature the muon relaxation is that expected for a sys- 
tem in the fast-fluctuation regime. The slow spin relax- 
ations observed in the frequency dependence of the ac 
susceptibility above 25 K are not within the muon time 
window so are not observed. The observation of two 
components in the muon relaxation below 24 K, together 
with a changing frequency dependence of the ac suscep- 
tibility, suggest a model of coalescing magnetic clusters 
forming well above Tjj. On cooling below ~ 50 K, these 
magnetic clusters occupy an increasing volume fraction, 
causing the increase in the frequency dependence of the 
ac susceptibility, until at ~ 25 K the sample possesses 



short-range antiferromagnetic order. Below 25 K these 
clusters coalesce, leading to the decrease in the frequency 
dependence of the ac susceptibility, and the increase in 
P s , until at Tn the sample possesses long-range order. 
Short-range order within the layers and antiferromagnet- 
ically correlated layers is consistent with the persistence 
of unbroadened magnetic Bragg peaks observed by neu- 
tron powder diffraction over a similar temperature range 
above Tf^Sa. 

In conclusion, NaNiC>2 shows the onset of long-range 
magnetic order at Tn = 19.5 K, with the dependence of 
the sublattice magnetization on temperature appropriate 
for a 2D XY magnet. The slowing of spin fluctuations 
above T s f is evident in the ac magnetic susceptibily data. 
At temperatures just above Tn there is evidence of short- 
range order, and of magnetic clusters persisting within a 
paramagnetic phase above this temperature. 
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